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Abstract 



We study flavour violation in the minimal SUSY SU(5) GUT assuming all the 
third generation Yukawa couplings to be due to the renormalizable physics 
above GUT scale. At large tan (3, as suggested by Yukawa unification in SU(5), 
sizable flavour violation in the left (right) slepton (down squark) sector is in- 
duced due to renormalization effects of down type Yukawa couplings between 
GUT and Planck scales in addition to the flavour violation in the right slep- 
ton sector. The new flavour physics contribution to K — K, B — B mixing 
is small but might be of phenomenological interest in the case of b — > 57. 
The sign of the latter contribution is the same as the sign of the dominant 
chargino contribution, thus making the constraints on SUSY scale coming 
from b — > S7 somewhat more restrictive. The most important feature of the 
considered scenario is the large rate of lepton flavour violation. Given the 
present experimental constraints, the \x — > ej and fj, — e conversion branching 
ratios are above the sensitivity of the planned experiments unless the SUSY 
scale is pushed above one TeV. 



I. INTRODUCTION 



The minimal supersymmetric standard model (MSSM) one of the best motivated 

extensions of the Standard Model (SM), has triggered intensive research both in theoretical 
as well as in experimental physics. Despite of non-trivial constraints on its parameter space 
coming from collider and low energy experiments the MSSM has so far successfully passed 
all the tests of precision physics. 

In the most general case MSSM contains more than one hundred free parameters. Some 
of them may give rise to unobserved phenomena like proton decay, large electron dipole 
moments, large flavour violation, etc. To explain the absence of such phenomena additional 
assumptions are needed to explain the pattern of supersymmetry (SUSY) breaking param- 
eters. An attempt towards that direction is to regard the MSSM as a low energy remnant 
of some grand unified theory (GUT) such as SU(5), which predicts unification of gauge 
couplings. Since the unification naturally occurs in the MSSM, one might try to apply a 
similar organizing principle to the soft SUSY breaking masses. Another prediction of GUTs 
is the unification of some or all Yukawa couplings of each generation. To achieve successful 
Yukawa unification [0-0] the ratio between the vacuum expectation values (vevs) of the two 
MSSM Higgs doublets, tan /?, is found to be in the range tan (3 ~ 30 — 50 for tau-bottom 
unification and tan/3 ~ 50 for tau-bottom-top unification ||. However, large mixings in the 
lepton sector may somewhat change this picture 0. 

Stringent tests of SUSY GUTs are offered by flavour violation experiments. The flavour 
violation in SUSY theories is not suppressed by the high scale where the SUSY breaking 
parameters are generated but rather by the mass scale of these terms themselves |§ which is 
believed to be of order TeV. In the scheme of the minimal flavour violation of Barbieri and 
Hall || the SUSY breaking parameters are generated above the GUT scale Mqut ~ 2 x 10 16 
GeV, at the reduced Planck scale Mp ~ 2.4 x 10 18 GeV by gravitational interactions and 
are therefore universal at Mp. The renormalization group (RG) evolution below Mp induces 
non-universalities of the soft terms at the GUT scale where the SU(5) gauge group breaks 
down to the usual MSSM gauge symmetry group. Therefore the flavour mixings present 
in the Yukawa couplings at GUT scale cannot be rotated away and should be reflected at 
low energies in the squark and slepton mass matrices. The phenomenology of this scenario 
in SU(5) SUSY GUT has been studied with the assumption that the top quark Yukawa 
coupling is the only sizable one |IJ] which implies flavour violation in the right-handed 



slepton and left-handed down squark sectors. Rates of the lepton flavour violating (LFV) 
processes fi —>■ e^y, fi — e conversion in nuclei and r — > yq are found to be large for some parts 
of the SUSY parameter space but due to the cancellations between gaugino and higgsino 
loops, they almost vanish for some other parts of the parameter space [O . 



Another set of SUSY theories predicting large rates of flavour violation are the ones 
with right-handed massive neutrinos [0,0- These models are motivated by the Super- 
Kamiokande results which imply maximal mixing between tau and muon neutrinos. The 
large mixing in the neutrino sector induces large flavour mixings in the left-handed slep- 
ton and right-handed down squark sectors in these models. The rates of flavour violating 
processes may in this case be much larger than in the minimal model, for example, the 
branching ratios of LFV processes can be close to or even exceed the present experimental 
bounds. 



2 



The aim of the present work is to revisit the minimal flavour violation scenario in the 
SUSY SU(5) grand unified theory. We extend the considerations of the previous papers 
T0 , I1J in several directions. Following the hints of possible Yukawa unification^ we consider 



the case when all third generation Yukawa couplings are large and given by the renormaliz- 
able physics above the GUT scale. In this case, sizable flavour violation in the left-handed 
(right-handed) slepton (down squark) sector is induced due to the renormalization effects 
of down type Yukawa couplings between GUT and Planck scales, in addition to the flavour 
violation in the right-handed slepton sector. Thus the pattern of flavour violation in SUSY 
SU(5) GUT at large tan/? resembles the flavour violation pattern in models with massive 
neutrinos, and deserves phenomenological studies. We do not make an attempt to identify 
and to study only that part of the soft terms parameter space in which tau-bottom Yukawa 
unification is achieved starting from the low energy parameters, but instead assume that 
there may be large corrections depending on the details of the GUT theory. We allow tan /3 
to be a free, although large, parameter so that the effects of r and b Yukawa couplings 
become non- negligible. The details of our calculations are given below. In our numerical 
analyses we consider the case where Yukawa unification is possible as well as the case where 
it is not. 

In addition, we calculate the new flavour physics contributions to K — K and B — B 
mixings as well as to b — > 37, fi — > ej, fi — e conversion in nuclei and r — > ^7 branching 
ratios. We also comment on the possibility of large GUT phases in this context. There 
is an increasing amount of constraints on the MSSM mass spectrum coming from direct 
collider searches as well as from indirect low energy measurements. We take these into 
account when calculating allowed ranges for the physical observables. In particular, the 
constraints coming from b — > S7 turn out to allow only the LFV processes to be in the 
phenomenologically interesting ranges. 

The main motivation for the present work is to study how sensitively the planned next 
generation LFV experiments will probe the SUSY model considered here. In the near future 
the branching ratio of the decay fi — > e'j will be probed with the sensitivity below 10~ 14 Jl5| 
and fi — e conversion in nuclei below 10~ 16 [JTB] (in more distant future \i factories may further 
reduce these numbers by orders of magnitude |T7[). We shall show that if these experiments 



will indeed reach the planned sensitivity then, in the scenario considered in this work, LFV 
must be observed unless the SUSY scale ^/m^m^ is above TeV scale. This conclusion 
remains valid even in the case where the only source of LFV is the right-handed slepton 
mass matrix as in Refs. [JTTJ since the deep cancellation in the branching ratios observed in 



Ref. [1 I] occurs at slepton masses which imply multi-TeV SUSY scale. In this context it 



is interesting to note that if we require approximate r-b Yukawa unification at GUT scale 
then the sparticle masses are required to be too high for direct production at future collider 
experiments. In this case, the observation of LFV may turn out to be the only signal of 
sup er sy mmet ry. 

The paper is organized as follows. In the second section we present the crucial parts 



1 Values of tan/3 between 0.7 and 1.8 are excluded by the direct LEP2 searches for the lightest 
MSSM Higgs boson @. At the end of LEP2 run, values of tan/3 below 2.6 will be probed Q. 
This constraint is not very restrictive but hints into the same direction as the Yukawa unification. 
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of model we are using throughtout this work, especially the relevant mixing matrices are 
introduced. In section III renormalization from the reduced Planck scale to the electroweak 
scale is discussed. In section IV we present our numerical results. Finally in section V we 
summarize our conclusions. 



II. SOURCES OF FLAVOUR VIOLATION 

In this section we present some details of the minimal flavour violation scenario considered 
in our work. We assume that the supersymmetric SU(5) grand unified theory is valid at mass 
scales between Mp and Mqut- I n the SU(5) model there are three generations of matter 
multiplets ipi and fa, % — 1, 2, 3, which form the 10 and 5* dimensional representations of 
SU(5), respectively, and 5 and 5* dimensional representations of Higgs multiplets, H 2 and 
Hi, respectively. The tenplets fa contain the quark doublets, the charged lepton singlets 
and the up-type quark singlets, while the down-type quark singlets and the lepton doublets 
are included in the fiveplets fa. The Higgs fiveplet Hi contains the MSSM Higgs multiplet 
Hi and a coloured Higgs multiplet H C i, and H 2 contains the second MSSM Higgs multiplet 
H 2 and another coloured Higgs multiplet Hci- An adjoint representation Higgs multiplet 
X causing the breaking of SU(5) should also belong to the Higgs sector of the model. We 
neglect the Yukawa coupling and the soft SUSY-breaking parameters associated with it. 
Inclusion of these terms will increase the non-universalities at GUT scale and thus increase 
the amount of flavour violation in the model. Thus the terms in superpotential W relevant 
for our consideration are given by 

W = \f u ^t B ^ D H^e ABCD E + V2f di .^t B fa A Hi B , (1) 

where A,B,... = 1,...,5 are the SU(5) indices. The relevant soft SUSY breaking terms 
associated with the SU(5) multiplets are 

- £SUSY breaking = {^V^iM^J + (m\)ij4>\<j>j + m^/lj^l + 777,^/4^2 

+ S-A Uij fafah 2 + V2A d ..fafah + h.c.\ , (2) 

where fa an d fa are the scalar components of the fa and fa chiral multiplets, respectively, 
and hi and h 2 are the Higgs multiplets. 

In the minimal SUGRA scenario the soft SUSY breaking parameters are generated by 
gravitational interactions and are universal at Mp: 

( m lo)ij = ( m i)ij = <% m o> 

.2 



™>hi = m h 2 = m , 



A-u^ — (A u • fu)ij, A diJ — (A d - fd 



'ij ■ 



A' Uij = A' di , = StjAo. (3) 



Below Mp till Mqut the parameters in Eq. (|I|) and Eq. (^) evolve with energy according to the 
renormalization group equations (RGE) of SUSY SU(5). These can be found for example in 
Ref. 0] and we do not present them here. Because the third generation Yukawa couplings 
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f U33 and fd 33 are much larger than the Yukawa couplings of the first two generations, large 
hierarchies in the parameters of Eq.(|3]) are induced at Mqut- Therefore the soft mass terms 
m^ , 777J and A' ud remain diagonal in the generation space but the third generation masses 
are smaller than the masses of the first two generations, which remain degenerate to a good 
approximation. 

To understand the origin of flavour violation in SUSY SU(5) GUT we have to discuss 
the evolution of the Yukawa coupling matrices fd and /„ with energy. At Mp we can choose 
the basis in which f u is diagonal, f p = (fu)a^ij- in that case the Yukawa matrix f d can be 
diagonalized with a bi-unitary transformation, 

w/y, (4) 

where V p corresponds to the Cabibbo-Kobayashi-Maskawa (CKM) matrix at Planck scale 
and P = diag(e i<?il , e 1 ^ 2 , e^ 3 ), where 01,2,3 are the GUT phases (satisfying 01 + 02 + 03 = 0) 
present in addition to the phase in V p . However, the mixing matrix U p can be rotated away 
by a redefinition of the fiveplet fields 0. No trace of that rotation will remain in the soft 
SUSY breaking terms because these are universal and proportional to unit matrix at Mp. 
Below Mp the Yukawa matrices fd and /„ run according to 

fuij + 6(fuf u fu)ij + 2(fdfdfu)ij + 2(fufdfd 



dt " 



lG7i 2 ^-fd.. 
dt 13 



-^ + 3Tr(/t/ u ) 
M gl + 4Tr(/]/ d ) f di . + 6(/jJ/ d ) y + 3(f u flf d ) tr (5) 
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Because the third generation Yukawa couplings are large, of order unity, the running of 
the off-diagonal elements of the Yukawa matrices is significant. Notice that despite of the 
chosen basis, non-zero off-diagonal elements are generated also in f u . To predict the running 
of the off-diagonal elements requires a precise knowledge of all the elements of the Yukawa 
matrices. This, however, goes beyond the assumptions made in our work: we assumed that 
only the third generation Yukawa couplings can be reliably estimated via the RGEs from the 
experimental data. Therefore, in our numerical estimates we assume that the off-diagonal 
elements of the rotation matrix U, which rotates the fiveplet fields are small, of order or 
smaller than the corresponding CKM matrix elements at GUT scale. 

At GUT scale the Yukawa coupling matrices f d and f u get renormalized to f® and 
respectively, and can be diagonalized again with bi-unitary transformations. However, now 
any rotation of the superfields if) and is reflected in the soft SUSY breaking parameters 
because these are not universal at GUT scale any more. These rotations give rise to large 
flavour changing effects, as will be discussed in the following. 

At Mqut the SU(5) gauge group breaks spontaneously into the usual MSSM gauge group 
SU(3) C x SU(2) L x U(l) y . The MSSM superpotential W valid below GUT scale is 

W = Q°{f Ui .)U°°H 2 - QHf^DfH, - Ef } {f e ..)L° i Hi - ^H X H 2 (6) 

and the soft SUSY breaking terms are given by 

- £ soft = Lfimi^L* + ^rK w + 4 0t KkQi + ur«ku? + TK)ii^ 

+m] i H\H l + m\ 2 E\E 2 + (Q°{A Uij )U? >H 2 - Q°(A d . j )Dj°Hi - Ef{A eij )L°Hi 
+B H H l H 2 + -M X BB + l -M 2 W a W a + ^M 3 g a g a + h.c)j . (7) 
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Here all the fields are explicitly written in the flavour eigenstate basis denoted by the su- 
perscript zero. The boundary condition for the soft terms are specified via 

= {ml)ij = (m|) y = (m? ) y , 

m \ Y = m h, , m H 2 = m h 2 , 

M x = M 2 = M 3 = M . (8) 



Since there is no splitting among the gaugino masses above the GUT scale, it is most 
convenient to stipulate M at the GUT scale. For the Yukawa couplings and trilinear terms 
this model predicts 

f? = f^ A? = A G (9) 

at GUT scale. 

From the GUT scale down to the electroweak scale the Yukawa couplings and the soft 
parameters evolve with energy via the MSSM RGEs [19]. Once the electroweak symmetry 
is broken, rotations of the superfields 

D° = V d D, E° = V e E, U° = V U U, 

D co = U* d D c , E co = U*E C , U co = U*U C , 1 ' 

bring quarks and leptons into their mass eigenstates with diagonal Yukawa couplings f' d , f' e 
and f' u 

f dl = {vjf d u:h, f ei = (utfeV e ) u , /; = (k t /^:)- (11) 

At the GUT scale the down quark and lepton masses are predicted to be equal. Therefore 
the diagonalizing matrices are related at Mqut as 

v d G = uf, u2 = v; G . (12) 

This implies that at low energies the left-handed quark rotation matrix, the CKM matrix in 
the basis in which up quark Yukawa matrix is diagonal, can be related to the right-handed 
lepton rotation matrix and the right-handed down quark rotation matrix can be related to 
the left-handed lepton rotations via the RGEs. 

At the same time the rotation in Eq.(|H]) changes the basis of the superpartners of quarks 
and leptons. For example, the mass eigenstates of the charged sleptons and sneutrinos can 
be expressed as 

where Te and T u are 6x6 and 3x3 rotation matrices, respectively. The slepton mass 
matrices are given by 

/ V^m L V e + m 2 e -m 2 z cos 2/3 (| - sin 2 9 W ) —fim e tan (3 + m e U%A'JU e 



ml 



\ -Ai*m e tan (3 + UlA' e U e m e \J\m 2 p \J e + m\ - m| cos 2f3 sin 2 9 W 



V}m\V e + im| cos 2/3)17,. (11) 
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and analogously for squarks. 

Notice that this definition of the diagonalizing matrices Te,d,u differs from the one 



originally given in Ref. [^D] where the matrices T relate the flavour eigenstates directly 
to the mass eigenstates. The advantage of the present notation in calculating the flavour 
violating observables in our scenario is the following. It allows us to perform the MSSM 
RGE running of the soft terms (which are diagonal in the flavour space) and the CKM 
matrix elements (which induce the flavour mixings) separately without constructing the 
squark mass matrices at high scales and without running each element of the 6x6 sparticle 
mass matrices. The slepton and squark mass matrices are then constructed at low energies 
in terms of the low energy values of the soft terms, quark masses and the mixing matrices. 

We use the standard notation for the neutralino and chargino mass matrices. The neu- 
tralino mass matrix in the basis (B, W°, H$, H$) can be diagonalized as 

M dia 9 = NM „ N i ( 15 ) 

where is a 4 x 4 unitary matrix. Similarly, the two diagonalizing matrices Ol,r in the 
chargino sector can be found from 

Mt 9 = L ( m M \ ^sin/W \ Q] 
x y v 2 cos pm w n J 

For the further details see, e.g. , Ref. M. 



III. RENORMALIZATION PROCEDURE 

Before calculating the rates of flavour violating observables, let us describe our pro- 
cedure of calculating the input SUSY parameters via the RGE evaluation. We start the 
RGE running at Mz where we introduce as the low energy input the values of the gauge 
coupling constants, the tau and bottom quark Yukawa couplings corresponding to the tau 
and bottom quark masses m T (M z ) = 1.784 GeV and m b (M z ) = 3.0 GeV [pl| , respectively. 
We parametrize the CKM matrix in the standard way |2^] with # 12 = 0.22, # 2 3 = 0.04 
and 6*i3 = 0.003. We evolve these quantities from Mz to mt using two loop SM RGE-s for 
five flavours. At m t we include top quark Yukawa coupling corresponding to the pole mass 
m t = 174 GeV and run the SM RGE-s for six flavours up to the scale Q where superparticles 
are introduced. At the scale Q we convert the SM couplings to the MSSM ones fixing the 
value of tan (3 and include the SUSY loop corrections to the bottom quark and tau lepton 
Yukawa couplings according to Ref. |4]j23[1 . Then the gauge and Yukawa couplings are eval- 
uated at GUT scale with the help of two loop MSSM RGE-s. For the running of the CKM 



matrix elements we use the one loop RGE-s from Ref. |24 . 

The crucial issue in the present context is the correct tau-bottom Yukawa unification 
as predicted by SU(5). The Yukawa unification has been studied extensively in literature 
HHJII and it turns out that unification is possible only for \i < and 30 ^ tan/3 ^ 50. This 
can be understood as follows. The correction to the bottom Yukawa coupling , 

A-T^r 1 (17) 

1 + b Vi 
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which is numerically the most important one, is proportional to tan (3 and its sign depends 
on sign(n). Tau-bottom Yukawa unification can be obtained if the top quark mass is ap- 
proximately at a fixed point at tan/? ps 1.6. On the other hand, small values of tan/3 are, 



anyway, excluded by the LEP 2 search for the lightest Higgs boson ||13yl4]| . To achieve 
the tau-bottom Yukawa unification therefore a sizable negative 5b is needed, which requires 
large values of tan/? and fixes the sign of the [i parameter to be sign(n) = — 1. On the other 
hand, it has been argued []5| that there might be important threshold corrections due to the 
coloured triplet Higgses at GUT scale modifying sizably the bottom quark Yukawa coupling 
already at Mqut- in any case, the coupling which receives large corrections is the bottom 
Yukawa coupling, tau Yukawa coupling is less sensitive to these corrections and thus the 
prediction of tau Yukawa coupling at GUT scale is more robust. Therefore, in our numerical 
calculations we use actually the matching condition 

/ d G [SU(5)]=/ T G [MSSM] (18) 

at Mqut to give a numerical value to the SU(5) Yukawa coupling fd 33 , and allow the value 
of the MSSM coupling f b G to be numerically different. To achieve semi-realistic Yukawa 
coupling evolution we fix the magnitude of 5b by hand for each considered value of tan/3 
but keep its sign to be determined by the sign(fj,). In such a case the Yukawa unification is 
achieved for sign(fi) = — 1 and not achieved for sign(fi) = +1. In our numerical examples 
we consider both possibilities. For definiteness we consider two values of tan/?, tan/3 = 35 
(corresponds to jf ~ 2/ G ) and tan/3 = 48 (corresponds to / G ~ / G ). 

At Mqut we fix the leptonic mixing matrices V e and U e according to Eq. flT2|) . Unlike 
the quark mixing matrices V e and U e do not run with energy due to the absence of neutrino 
Yukawa couplings in the MSSM. Choosing the basis in which the top Yukawa matrix is 
diagonal we have Vd = Vckm and therefore \U e \ij = \V G KM \ij. There is no experimental 
restriction on the mixing matrix of the right-handed quark fields. Therefore we assume that 
large couplings fd 33 , f U33 in Eq.(|) generate small non-zero (31) and (32) elements for Ud and 
V e but the angle #12 remains negligible in these matrices. Numerically we consider three 
cases: = U^ 32 \ = 0.1 x C/( 31 )-( 32 ) and = 0. 

Further, we run the SU(5) gauge and Yukawa couplings from Mqut to Mp where we fix 
the values of mo and Aq. For the trilinear coupling we take always Aq = in our numerical 
calculations. Then we run all the parameters, including the SU(5) soft terms down to the 
GUT scale, apply the boundary conditions Eq.(g) and run all the MSSM parameters down 
to the scale Q which we take to be 200 GeV. At that scale we fix the SUSY parameters 
[i and B via the conditions of spontaneous symmetry breaking and calculate the masses 
and mixing matrices of the sparticles. The chargino and neutralino mass matrices are given 
by the standard expressions which can be found for example in Ref. M. In order not to 
run into contradiction with the mass bounds on the SUSY particles from direct searches 
we take M <: 150 GeV and m ^ 150 GeV; the lightest sparticle masses corresponding 
to m = M = 150 GeV and tan (3 = 48 are roughly M fl = 73 GeV, M-+ = 98 GeV and 
M^o = 58 GeV and are just on the limit of the present LEP bounds. 



IV. RATES OF FLAVOUR CHANGING PROCESSES 
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A. Kl — Ks System 



We start our numerical analyses by calculating the new SUSY flavour changing physics 
contribution to physical observables in the K^ — Ks system. Recently it has been emphasized 



25fl that in models with massive neutrinos the SUSY flavour changing contribution to €k 
may be large, especially if new GUT phases are present which may maximize the effect. 
Because at large tan/? the pattern of flavour violation in the SUSY SU(5) is the same as in 
the models with right-handed neutrinos we study this issue carefully. 

The dominant new physics contribution to AS = 2 processes comes from box diagrams 
with gluinos and down squarks running in the loop |20|| . Because in our scenario the first two 
generation sfermions are almost exactly degenerate the appropriate tool to use is the mass 
insertion approximation | ]2G|| . The current state of the art on this subject is summarized in 
Ref. |27j in which the low energy AS* = 2 effective Hamiltonian is calculated including NLO 
QCD corrections and the relevant hadronic matrix elements are evaluated using the lattice 
results for the -B^-parameters. We use the expressions for the AS* = 2 Wilson coefficients, 
NLO QCD corrections and hadronic matrix elements as well as numerical input exactly as 
in Ref. p7|] and therefore we do not copy them here. However, we need to specify the model 
dependent input. 

In the mass insertion approximation the flavour violation is characterized with the di- 
mensionless parameters (Sll)^, (Shr)^, (^LR)ij and (Srl)^ defined as 




J_ ( {ml) LR 
fh 2 \ {m%)RL (m%)RR 



(19) 



where (m 2 ^)MN, M, N = L, R are the corresponding 3x3 components of the 6x6 down 
squark mass matrix and fh 2 is an average squark mass appropriate for the problem under 
consideration. We calculate the down squark mass matrix as described in detail in previous 
sections. The mixings among the left down squark fields are generated by the CKM matrix 
Vckm- The mixings among the right down squarks are assumed to be given by the matrix 
Ufkj — (YcKiv^iji i->j = 2,3 with vanishing angle 6*12. Thus the off-diagonal elements of 
are smaller than the elements of Vckm- Because only the (12) components of 5mn enter to 
the AS = 2 process we identify fh 2 = (m%)n. 

The Kl — Ks mass difference and the CP-violating parameter ex are given by: 

AM K = 2 Re(K°\H* s f = 2 \K } , 

- = VI^ ^ K °^r\K°) - (20) 

It turns out that the new physics contribution to AMk is very small in our model, typically 
at the permil level of the measured value. However, the new contribution to ex which is a 
small quantity in the SM can be sizable. In Fig. |we plot the gluino mediated contribution 
to €k as a function of the average squark mass fh fixing the gluino mass to be Mg = 420 
GeV and assuming that the new GUT phases maximize the effect. For the squark masses 
of order 500 GeV and large tan (3 the contribution may exceed 25% of the measured value 
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FIG. 1. Gluino mediated contribution to ex as a- function of the average down squark mass 
assuming that the possible new GUT phases maximize the effect. 



ex = 2.3 ■ 10~ 3 . However, because of the assumption of the GUT phases and because of 
large theoretical errors in the hadronic matrix elements no useful constraints can be derived 
from €k on the parameters of the model. The reason why the contribution to ex in SUSY 
SU(5) is smaller than in the models with massive neutrinos is the smallness of the (5 mn) 12 
elements. The magnitude of the largest mixing parameter (£ll)i2 is typically below 1CT 4 
and for the chosen parameters (5rr) i 2 is about factor of three smaller. The (LR) parameters 
are smaller by three orders of magnitude. We conclude that the small mixing angles in our 
approach suppress the contribution to the K — K mixing. 



B. B d - B d System 



Here we estimate the new physics contribution to the B d meson system. The SM con- 
tribution to AM Bd = 2\M 12 \ is given by |2§ 



M 



SM 
12 



12^2 



VQCDB Bd fiM Bd M^(V td V t tyS (z t ) 



where 



S {z t ) 



Az t -llzl + z\ 

4(1 - z t y 



"iz\ In z t 

2(1 - z t y 



(21) 



(22) 



where z t = m 2 t /m 2 w andr/ QCD = 0.55±0.01. With B Bd = 1.29±0.08±0.06 and f Bd = 175±25 
MeV it successfully predicts the measured value AM Bd = 0.470±0.019 ps -1 . To estimate the 
magnitude of the dominant gluino induced contribution to B — B mixing in our scenario we 
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use again the mass insertion approximation. Neglecting the small LR and RL contributions 
it reads [EST] 



M 



SUSY 
12 



o; 



216m 2 3 



l -B Bd f B M Bd {((5 LL ) 2 31 + (5 RR ) 2 31 )(mf 6 (x) + 2Axf 6 (x))+ (23) 



(S L Lhi(S RR ) 3 i ^36-24 
= 67^ ( 



B d 



\m b + m d 



f 6 {x)+ 72 + 384 



/ M 



\m b + m d 



xf&{x) 



x 3 - 9x 2 - 9x + 17 + 6(1 + 3z)lnx 



3(x 



-x 6 - 9x + 9x + 1 + 6x(l + a;) hue) . 



(24) 
(25) 



where x = M~/rh . In this case the common squark mass is taken to be in = (mj~)33 and 
the mixing elements £31 are calculated as in the previous subsection. 

In the case where the sparticle masses are assumed to take their lower allowed limits, the 
SUSY contribution to AMs d turns out to be at most of order a few percent. Taking into 
account the large errors in the input parameter Jb Bs d no useful constraints on our model 
can be derived drom the measurement of AMb,- 



The radiative decay b —>■ 57 is known to get large contributions from SUSY particle 
loops and therefore this process implies strong constraints on the allowed SUSY parameter 
space pH|f2D|j3l]|] . Besides the SM PU-boson-i-quark contribution there are also charged Higgs, 



chargino, neutralino and gluino contributions. At the electroweak scale all of them have been 
calculated in Ref. PD| . While in the SM the NLO QCD corrections are included to b — * 57, 
in SUSY theories the NLO analyses has been performed only for specific scenarios [^TJ. The 
complication here lies in the new flavour structure of SUSY theories, as compared to the SM, 
which extends the operator basis beyond the SM one. Therefore the LO QCD corrections 
to the gluino contribution which gives the dominant new flavour physics contribution was 



calculated very recently in Ref. [32|. Before that, the gluino contribution in a generic SUSY 



flavour model was studied in [26| without including QCD corrections. 



Here we shortly review the results of Ref. The effective Hamiltonian for b — > 37 is 
expressed in two parts 

Keff = H e ff M + ft 9 effi ( 26 ) 

where T~leff M ^ s ^ e effective Hamiltonian in which the structure of flavour violation is the 
same as in the SM and the gluino contribution 7i 9 e ^ exhibits the new flavour structure. The 
Wilson coefficients of the first term, 

AC 1 

n c e « M = —^V tb V* £ C^)0^) , (27) 

contain the SM as well as the charged Higgs, chargino and neutralino contributions. The 
dominant operators in Eq. (|27j) are the dimension six magnetic operators (fn b (fi) is the run- 
ning mass) 
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° 7 = TEr? F ^ ' ° 8 = ife mb ^ (sa^T a P R b) G% , (28) 

and the operators 0' 7 8 obtained by L <-> i?. The gluino effective hamiltonian 

= E <3»o<» + EE <WoWaO • (29) 

i i q 

contains in addition to the dimension six magnetic operators On,g, ^n,gi ®8b,g, &sbg a ^ so 
operators of dimension five, 

Org, 9 = eg 2 s (fi) (sa^P R b) F, v , 8M = g s {n) g 2 s {n) (sa^T a P R b) G% , (30) 

in which the chirality-violating parameter is the gluino mass in the corresponding Wilson 
coefficients, as well as new four-quark operators Of g , q — u, d, c, s, b, which are listed in 
Ref. [32|. These new operators change the structure of the LO QCD corrections in general 
SUSY flavour models. 

The Wilson coefficients in Eq.fl2"7|) including the LO QCD corrections in the MSSM are 
well known. Their explicit expressions can be found in |30] and we do not rewrite them here. 
Instead, let us for a moment concentrate on the study of the effective hamiltonian Eq.([29"D 
in our model. 

At the electroweak scale the relevant Wilson coefficients are given by |32] 



6rf 6 1 

C 7b ~ g {nw) = -yrA C(R) 53 — T { T dl t *dl) /: 

k=l d k 



2\ x gd 



k) ' 



Crgdnw) = M 9 -^- 2 c(R) £ ±r (r% R r* D k £) f. 

fc=l d k 



&{%gd t 



k) 1 



<WM = -T^E-V(r^r^) [{C{R)-±C{G))f 2 {x gdk )-\C{G)h{x gdh )\ , 
ion k=1 

Cmfaw) = m 9 — 53 _L (r» fl r;&) [(C(r) - jc(go) um - \c{G)h{x 9dk )\ , (31) 

k=l d k 

Here the matrices Te>l and are the 6x3 submatrices of T^, 

F 6x6 _ ( -n6x3 p6x3 \ /oo\ 
L D y DL DR ' J [ 6Z ) 

and the ratios x g d k are defined as x g d k = M~/m~ . The Casimir factors C(i2) and C(G) are 
C(i?) =4/3 and C(G) = 3 and the functions fi(x), i = 1, ...,4, are given by 

fi(x) = — — (x 3 — 6x 2 + 3x + 2 + 6xlogx 



12 (x - 1 



M*) = 7777 1 1U (2^ 3 + 3x 2 - 6a; + 1 - 6x 2 log 
12 (x — l) 4 v 

/s(x) = 3 (a; 2 - 4s + 3 + 2 logo; 

/ 4 (x) = 1 (a; 2 - 1 - 2a; log a;) . (33) 
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The Wilson coefficients of the corresponding primed operators are obtained through the 
interchange «-> T^ L . At the low scale /j, b , the LO renormalization group improved 

coefficients become 



27 8 / 25 27 \ 

C 7 b,§{Hb) = rfi C 7bi g(n w ) + | (77S - r^i) C 86 , § (// w ) + R 7b ,g(/i b ) 



(34) 



where r/ = a s (yUvF)/o ; s(A i fe)- The remainder function R 7 b, g (Hb) turns out to be numerically 
very small and we neglect it in our numerical computation. The low-scale Wilson 
coefficients for the corresponding primed operators are obtained by replacing in (|34"|) all the 
unprimed coefficients with primed ones. 



A C,gu 



A C,gMn 



M,=420 GeV 
tan 0=35 



sign fi= + 



sign n= — 



800 1 000 

mJGeV] 



U °- 15 
< 



M g =420 GeV 
tan £=48 



sign /[/= + 



sign \jl=- 



800 

m sh [GeV] 



FIG. 2. Gluino mediated contribution to the Wilson coefficient C 7 as a function of the sbottom 
mass for fixed gluino mass M g = 420 GeV and tan/3 = 35, 1^8 as indicated in figures. The sign of 
the contribution depends on sign(p). 



The decay width of b — > 57 can be written 

m 5 b G 2 F \V tb V*\ 2 a 



where in our model 



Here 



C 7 9 



T(b -> try) 



32tt 4 



C e 7 ff 



C 7 



eff 



c 7 + c 9 7 



+ 



c 7 9 



c 9 7 = - 



16V2ir 3 a s (ii b ) 
Gf V t bV t * 

lQy/2Tc 3 a s (ij, b ) 
Gf VtbVts 



1 



C 7 b,g{^b) H C 7 g ; g(fM b ) 

m b 

C 7 b,g{^b) + —C 7 g-g{^b) 

m b 



(35) 
(36) 



(37) 
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and C 7 stands for the total contribution from the effective Hamiltonian Eg. (p7[) . In the SM 
its value is C 7 5M = -0.3 @. In the MSSM C 7 receives large contribution from charged Higgs 



loops which add constructively with Cf M and at large tan (3 also from chargino loops which 
add constructively (destructively) with Cj M for sign(fi) = — 1 (sign(fi) = +1). Therefore 
for large tan f3 and small Higgs and sparticle masses the C 7 may be completely dominated 
by SUSY. Therefore the present experimental result 

2 x KT 4 < BR(5 -> X s7 ) < 4.5 x 10" 4 , (38) 

which favours the following allowed range: 

0.25 < |C 7 e// | < 0.375, (39) 

implies strong constraints on the SUSY masses. 

Let us now study the gluino contribution to C 7 numerically. In Fig. ||] we plot the 
values of C 7 as a function of the lightest bottom squark mass m~ b for a fixed value of the 
gluino mass M g = 420 GeV and for two values of tan/3 and sign(fi), as indicated in the 
figure. Note that this is the smallest gluino mass consistent with chargino mass bounds 
in our model. The values of Cf might be sizable for small sparticle masses but decreases 
rapidly if the masses are higher. The important behaviour to notice is that the sign of C 7 
depends on the sign(fi) exactly the same way as the dominant contribution C-j. This implies 
that C7 and C 7 add up constructively and no cancellation between them is possible. It has 
been argued in Ref. [|3~2"| , |3~4"f that the constraints on the SUSY parameter space coming from 
C 7 which are rather restrictive can be relaxed by new flavour contributions. However, in our 
model this does not happen and the h —>■ 37 bounds cannot be relaxed with help of C 7 . 

The gluino contribution Cf is induced by the mixings in the left-handed down squark 
sector, thus by Vd = Vqkm- However, the coefficient C/ is almost entirely induced by the 
mixings in the right-handed down squark sector. If these mixings are similar in size to 
the CKM mixings, as we argue here, then in our scenario with large tan/3, C 7 is as large 
as Cf. While its absolute contribution to C 7 is subdominant (see Fig. 0) it still may 
have important phenomenological consequences. Namely, while the mixing induced time 
dependent CP asymmetry in B — > M s j, where M s is some CP eigenstate, is predicted to be 
very small, a few percent, in the SM; our model, where it can be expressed as 

T(t) - T(t) 2 Im \e- ieB i (C 7 + C$)C'f 
Lj! — A t sin AM Bd t , A = L -4J- 7> ?J , (40) 



r(t)+f(t) d ' |c 7 e// | 2 

with Q^i — a r 9{Mi2 d ) being the phase in the B — B mixing amplitude, asymmetries A t of 
more than 10% (which would be a clear and powerful signal of beyond the Sm physics) are 
allowed. Similar conclusions hold also in models with right-handed neutrinos fl35|. 



D. Lepton Flavour Violation 

So far we have shown that in SUSY SU(5) the new flavour physics contribution to flavour 
changing hadronic observables is subdominant. At the same time the SUSY contribution to 
b — ► 57 induced by the effective Hamiltonian Eq. (|27|) constrains severely the SUSY scale in 
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our model. Now we turn to study the LFV processes which are completely dominated by 
the new physics. The amplitude for the process lj — > Zj 7 * where the photon is off-shell can 
be written as 



M = ee a % 



q 2 la (A[P L + A R P R ) + m h ia aP q^ [A L 2 P L + A R P R )\ I, , (41) 



where q is the photon momentum and A\' R are the form factors giving rise to the process. 
Note that all the form factors contribute to fi — e conversion but only A 2 ' r give rise to the 
lepton radiative decays. The form factors are induced by two types of loop diagrams with 
neutralinos and charged sleptons in the loop, and charginos and sneutrinos in the loop: 

A[f = A ( $ L > R + A$ L ' R , (42) 

where schematically written 

A^ L ' R = A^ L ' R (M^,m h N,T E ) , 

4f' R = 4f"* , mp, O l , r , r„) , (43) 
depend on the slepton masses and mixings as well as on the neutralino and chargino masses 



and mixings. We have adopted the formulas for the form factors from Ref. |L2| and we do 
not present them here. The decay rate of lj — > k'j is then given by 



r(/^/ l7 ) = ^-mf 3 (|^| 2 + |^| 2 ) . (44) 

We start with studying the decay \x — > e 7 in the case where the only source of LFV is 
the mixing in the right-handed slepton sector as given by Eq.fll^), and the mixing matrix 
V e is the unit matrix. This is the case studied in Ref. |T| . In Fig. [| we plot the branching 
ratio of fi — > cy on the plane of the lightest charged slepton mass mi and the lightest 
chargino mass M$+ for two values of tan/5 = 35, 48 and sign(n) = —1 (to achieve Yukawa 
unification). We have taken into account the experimental constraints coming from b —> sj 
by requiring that the total value of C% is in the allowed range (551) . As seen in Fig. 131 the 



y 7 

constraint Eq. (j3~9D puts strong lower bounds on the lightest sparticle masses. Our values of 
the branching ratio of /i — > e 7 are much smaller than the quoted values in Ref. The 
reason is twofold. First, the b — > s 7 constraint pushes the sparticle masses to high values^ 
and this has not been taken into account in Ref. [11]. Second, the authors of Ref. fnj fix 



the top Yukawa coupling at Mqut to be / t =1.4 which for large tan/3 implies by far a too 
large top quark mass. For m< = 174 GeV and tan/3 = 35, the correct value is / t = 0.56. 
For large part of the parameter space the destructive interference between the gaugino and 
higgsino contributions suppresses the \x — > e 7 branching ratio to almost vanishing values. 



2 It has been noticed in Ref. |36| that for a particular corner of the parameter space where Mo -C 
rriQ ~ A one can supress b — > s 7 and still have light gauginos. However, also LFV processes are 



suppressed for this parameter space. 
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\i — > e y 



|X->ey 




1300 



FIG. 3. The branching ratio of fi — ► e"f as a function of the lightest charged slepton and the 
lightest chargino mass assuming that the flavour mixing in the left-slepton sector is vanishing, 
V e = 1. We have fixed sign(p) = — and values of tan/3 are indicated in the figures. The destructive 
interference between the gaugino and higgsino contributions suppresses the branching ratio in large 
regions of the parameter space. The constraints coming from b — > S7 are taken into account. 

The situation changes completely if we allow also flavour mixings in the left-slepton 
sector and allow sign(fi) to be also positive. Taking = U l J = {VQ Khl )^ as we predict in 
our model we plot the /1 — > branching ratio for tan j3 = 35 in Fig. [| and for tan (3 = 48 
in Fig. |5|. The branching ratios are about two orders of magnitude higher than in Fig. |3| 
and no cancellation occurs for any sparticle masses. This implies that fi — > ej is dominated 
by the sneutrino-chargino contribution. The LFV pattern here is exactly the same as in 
models with right-handed neutrinos. Notice that for sign(fi) = +1 the sparticle masses are 
allowed to be much smaller than in the other case. This is because for sign(fi) = +1 the 
chargino contribution to b — > 57 interferes destructively with the SM and charged Higgs 
contributions. In particular, for tan (3 = 48 and for very small chargino and slepton masses 
the chargino contribution can be so large that it cancels the SM and H + contributions and 
gives the allowed C7 value with an opposite sign. This is seen in Fig. |5| for sign(^) = +1 
in which a small parameter region around M^.+ « 100 GeV and w 300 GeV is allowed. 
This region can be excluded by collider searches for a very light chargino or by improving 
the bound on the \x — > ej branching ratio by a factor of few. 
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\i — > e y 




1 100 



FIG. 4. T/ie branching ratio of fi — ► e7 as a function of the lightest charged slepton and the 
lightest chargino mass for fixed tan/3 = 35. XTie sign of the u "parameter is indicated in the figures. 
Here and in the following = = (Yckm)^ ■ The constraints on SUSY mass spectrum coming 
from b — > S7 are taken into account. 



(1 — > e y |i — > e y 




1200 



FIG. 5. T/ie same as in Fig. |^ 6ut /or tan/3 = 48. 
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We have shown that the non-vanishing flavour mixings in the left-slepton sector are 
crucial to ensure the detectibility of fi — > ey in the planned experiments. Since there is 
no experimental information on these mixings the central question to ask now is how small 
the off-diagonal elements of V e can be and still allow successful determination of // — > ej. 
To analyze this question we plot in Fig. |^ the fi — > c-f branching ratio against the lightest 
slepton mass mi for a fixed M2 = 460 GeV which is roughly the minimal chargino mass 
for sign(fi) — — 1. The curves denoted by a correspond to V e = 1, curves denoted by b to 
= 0.1 x U % J , % 7^ j and curves denoted by c to = U l J . As can be seen, if the off 
diagonal elements of V e are as small as 10% of the corresponding U e elements then the deep 
cancellation is superseeded. Let us also mention that for the chosen chargino mass in Fig. 
P the SUSY scale, Msusy = y/ m ii m i 2 i * s ^susy ~ 1200 GeV for the minimally allowed 
slepton mass. Thus the LFV processes are sensitive to SUSY scale above TeV. 
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FIG. 6. The branching ratio of \x — > e^f as a function of the lightest charged slepton mass for 
fixed wino mass M2 = 460 GeV. We have fixed sign(/j.) = — and values of tan (3 are indicated in 
the figures. For the curves denoted by (a) V* 3 = 0, i 7^ j, thus they correspond to Fig. [|. For the 
curves (b) = 0.1 x U l J and for the curves (c) = U l J . 



Our calculations show that the rate of /1 — e conversion in nuclei is about 6 x 10 -3 times 
the branching ratio of /i — > cy. The qualitative behaviour of the /i — e conversion rate with 
the sparticle masses is the same as in the case of /i — ► cy. Thus the results for fi — e conversion 
can be obtained by rescaling the figures for /i — > cy. Therefore we do not present any new 
plots for the /1 — e conversion process. We conclude that the planned fi — e conversion 
experiments [^] are as sensitive to our models as the planned fi — > ej experiments ||15| . 



Finally let us discuss the decay r — > /ry. In Fig. ^ we plot the branching ratio of the decay 
r — ► /ry for the same choice of parameters as in Fig. ^ Even for tan (3 = 48 the branching 
ratio is always below a few times 10 -9 and unobservable in the planned experiments. Thus 
if r —>■ fi'-f will be discovered at these experiments this implies some other LFV scenario than 
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the one considered in this work. 




400 600 800 1000 1200 400 600 800 1000 



mJGeV] mJGeV] 
FIG. 7. The same as in Fig. but for the decay t —> fij. 
V. CONCLUSIONS 

Motivated by the sensitivity of the running or approved experiments to flavour violating 
processes we have studied flavour violation in the minimal SUSY SU(5) GUT at large tan (3. 
In this case the flavour mixing occurs both in the left and right slepton and squark mass 
matrices and are enhanced by the large value of tan (5. We have calculated the new physics 
contributions to K — K and B — B mixings and to the decays b — > 57, /x — > ej, r — > /ry 
and to fi — e conversion in nuclei. To predict reliably the rates of these processes we have 
correctly taken into account the measured values of the low energy parameters as well as 
the constraints on the SUSY particle masses. 

We found that in our model the new physics contributions to AMk and AM# are 
negligible, but might reach a 10% level in ex if there exist new GUT phases. No useful 
constraints on the model parameters can be derived from these processes. 

The decay b — *■ S7 receives contributions from two sources of flavour violation: from the 
loops proportional to the CKM matrix elements and from the loops exhibiting new flavour 
violation in the squark mass matrices. The latter contribution interferes constructively with 
the dominant chargino contribution. At large tan j3 the experimental constraints on the b — > 
57 branching ratio imply stringent constraints on the SUSY particle masses, especially for 
sign(fi) = — 1 as required by Yukawa unification. In this case, the SUSY scale is constrained 
to be at least TeV. For such a high squark masses the new flavour physics contribution to 
b — > 57 branching ratio is a few percent. Nevertheless, this may induce CP asymmetries 
considerably larger than in the SM. 

There is a competition between the sensitivity of the future LFV experiments to the new 
flavour physics and the constraints on the SUSY scale coming from the b —>■ 37 branching 
ratio. If the branching ratio of the decay \i — > cy will be tested down to 10~ 14 and the SUSY 
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scale is below 1 TeV then, the present scenario predicts that /i — > &-/ should be discovered 
in these experiments. The branching ratio of the decay r — > /ry is, however, predicted to be 
below a few times 1CT 9 and should not be seen at LHC in the minimal SUSY SU(5). 
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